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Abstract. - We study the onset of the wave-resistance due to the generation of capillary- 
gravity waves by a partially immersed moving object in the case where the object is hold at 
a fixed immersion depth. We show that, in this case, the wave resistance varies continuously 
with the velocity, in qualitative accordance with recent experiments by Burghelea et al. (Phys. 
Rev. Lett. 86, 2557 (2001)). 



Introduction. - The dispersive properties of capillary-gravity waves are responsible for 
the complicated wave pattern generated at the free surface of a still liquid by a disturbance 
moving with a velocity V greater than the minimum phase speed Vc — (4g7/p)^^^, where g is 
the gravity, 7 is the surface tension and p the density of the fluid |0] . The disturbance may 
be produced by a small object partially immersed in the liquid or by the application of an 
external surface pressure distribution The waves generated by the moving perturbation 
propagate momentum to infinity and, consequently, the disturbance experiences a drag R 
called the wave resistance [||. For V <Vc the wave resistance is equal to zero since, in this 
case, no propagating long-range waves are generated by the disturbance [Q . 

A few years ago, it was predicted that the wave resistance corresponding to a surface 
pressure distribution symmetrical about a point should be discontinuous a.tV — Vc § . More 
precisely, if Fq is the the total vertical force exerted on the fluid surface, the wave resistance 
is expected to reach a finite value i?c > for V . For an object much smaller than the 

capillary length = \/ ^ j pg, the discontinuity i?c is given by: 

i?c=-^-. (1) 

Experimentally, the onset of the wave-resistance due to the generation of capillary-gravity 
waves by a partially immersed moving object was studied recently by two independent groups 
While Browaeys et al. ||] indeed find a discontinuous behaviour of the wave-resistance 
at y = Vc in agreement with the theoretical predictions, Burghelea et al. observe a smooth 
transition. 

The discrepancy between the theoretical analysis of Q and the experimental results of 
might be due to the fact that the experimental setup of Burghelea et al. uses a feedback 
loop to keep the object at a constant depth while the analysis of |^ assumes that the vertical 
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Fig. 1 - We study the behavior of a hydrophobic moving object of characteristic size b immersed at a 
depth h. The velocity of the object is V. We decompose the force exerted by the fluid on the object 
in a component orthogonal to the free surface (Fo) and a component parallel to it (R). 



component Fq of the force exerted by the disturbance on the fluid does not depend on the 
velocity V (we might call such an analysis a "fixed force" analysis). In order to check this 
proposition, we perform in this letter a "fixed-depth" calculation of the wave-drag close to 
the onset threshold. A somewhat similar analysis was performed in the large velocity limit 
by Sun and Keller We will show that such a calculation indeed yields a cancellation of 
the vertical force Sut V = Vc, that is to say, according to equation (|l|), a smoothing of the 
discontinuity. 

Model. - We take the x, y plane as the equilibrium surface of the fluid. The immersed 
object exerts a stress at the fluid surface that can be considered equivalent to a pressure field 
p ^ that travels over the surface with a velocity V in the x direction. We assume p (the 
Fourier transform of p) to be of the form: 

)=i^o0(fc), (2) 

where k = \Jk'^ + ky and (^(0) = 1. In this case, p is isotropic and Fq is the total vertical 

force exerted on the fluid. 

Within the framework of Rayleigh's linearized theory of capillary-gravity waves, the Fourier 
transform ^(k) of the free surface displacement f (r) is related to the pressure held through pT|: 



k I (f){k) 



where Wo(^) = 5^ + l^^/p is the free dispersion relation, = fc^ — ik • V/v and v = rj/p is 
the kinematic viscosity of the fluid. 

Let us suppose the object is located at the origin of the moving frame. If h is its depth, 
the free surface displacement ^ must he —h under the pinpoint (here we suppose the object 
sufficiently hydrophobic and h not to large so that the pinpoint does not pierce the surface). 
This leads to the following normalization condition: 

The value of /i as a function of Fq can be readily calculated using equations (|^) and (||): 

h = E{V)Fo, (5) 



with 
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EiV) ^ f^^J m ] (6) 

Finally, the drag-force R is calculated by simply integrating the pressure force over the 
free surface [Q. 

R = - 1 d2rp(r)Ve(r) = - j ^ ^kp* (k)|(k). (7) 
This yields, using the explicit expression (||) for ^, 

R = i^o'A(F), (8) 

with 

MV) ^ ^ (9) 

According to equation (||), the integral A describes the fixed- force behaviour of the wave- 
resistance. Due to the symmetry of A is parallel to V and we shall henceforth set A — Au, 
where u — Y /V is the unit vector parallel to the velocity of the object. The authors of 
studied the properties of A in the case of a non-viscous fluid for which they showed that: 

(a) A = for V < V^; 

(b) A is discontinuous for V — > , with 

lim A = Ac = ^-; (10) 
v^v+ 2V27 

(c) in the large velocity limit 

In the case of a fixed depth analysis, Fq becomes a function of V . Using equations (||), we 
can rewrite the wave-resistance as 

In general we have to rely on numerics to calculate the integrals S and A. Typical results 
are presented on Figure |^ for an objet of size 0.1 mm immersed in water, and for a step-like 
function equal to 1 for k < 1/b and otherwise. 

We first observe on Fig. (||-a) that S increases sharply near the threshold. This leads 
to two rather different behaviours for A and R as shown on Fig. (^-b): while A exhibits a 
discontinuity close to V = Vc |jl^, the fixed-depth wave-drag R — K/'E? cancels smoothly at 
the critical velocity. 
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V/Vc V/V^ 

Fig. 2 - (a) Numerical calculation of the dimensionless integral H — H('K)/H(0). (b) Comparison of 
K{V) = A(V)/Ac (dashed line) and R — A/H^ (full line) describing respectively the fixed-force and 
fixed-depth behaviour of the wave- resistance. We clearly see that contrary to the fixed-force analysis, 
the fixed-depth calculation does not yield any discontinuity at the threshold. Here u = 10~® m^ • s~^, 
7 = 72 mN ■ m~\ p = 1000 kg ■ m^^* and b = 0.1 mm. 



Inviscid flow. - The characteristic features displayed by the plots of Fig. ^ can be 
captured by a zero-viscosity analysis. Setting i/ = 0, equation (|^) can be simplified as: 

rj^fc/ m ^ (13) 

where V denotes the Cauchy principal value of the integral. The integral is calculated in polar 
coordinates {k,6), where 9 is the angle of k with respect to V. Introducing the function G 
defined by: 

d6l / 1 



where mf. = k/ k + n/k and M = V/Vc is the "Mach" number, equation ( [l3| ) can then be 
rewritten as: 

1 dk ~ 

E^-l -mCik). (15) 
Using the residue theorem [Q, we get: 

where 9 the is the Heaviside step-function. The variations of with k are plotted on Fig. ^: 
reaches its minimum value (TO|)min = 2 for k = k. It shows that equation 2M^ — m| = 

has two solutions ki and ^2, with fci < k < fc2, if is larger than the critical velocity Vc, and 

none if V < V^. 

For V > Vc, 'E. evaluates to 
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Fig. 3 - Variations of m| (see text) with k. Since m| reaches a minimum for k = k with = 2, 
equation m| = 2M^ presents solutions for > 1. In this case, there are two possible values fci and 
k2 with fci < ft < ^2. 



1 



dk 



m 



0(fc) 



(17) 



The above integrals can be calculated in the two limiting cases A4 ^ 1 {i.e. V Vc) and 
M > 1 {i.e. F> t4). 

For large A4, the integrals are restricted to either large or small values of fc. If vanishes 
faster than 1/k for large k, we can show that the small k contribution dominates. In this 
region, the dispersion relation is dominated by gravity waves, so that we can approximate rri^ 
by k/K. A straightforward calculation then yields the following asymptotic expansion for S: 



1 

6717 



(18) 



Combining this result with (hWj and [13, we get 
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R - Unph- 



V 



Let's now focus on the case V close to Vc. We set: 



ki 



dk 



m 



2tt ^mi{mi-2M^y 



and 



dfc 



Ik, 27r ^mi{mi-2M^) 
Using the fact that ki and k2 are the roots of — 2M^, we can rewrite Si as: 
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■=1 



fci 



dfc k(P{k) 



"^TT ^{kyK^){k-kl){k-k2) 



When V Vc, the (fc — fci)(fc — fc2) term cancels in fc = k. Since all other terms are 
regular, we can write at the leading order: 



K- 



<P{k) dfc 1 



V2 Jq 27r ^(fc-fci)(fc-fc2)" 
This latter integral is readily calculated and gets a very simple form in the limit V Vc- 

S - In (fc2-fci)--^ In 



V2 ' ' 2V2 

Since is invariant by the transformation k/n — i- /t/fc, we see that Si and S2 have the 
same asymptotic behaviour for V Vc- In this limit we have S = + ^2^ /7K ~ 2S1/7K, 
hence: 

E^^^(.)ln(A.-l). (19) 

In the case of an object of size b, the width of (p is about 1/6. If we choose b much smaller 
than the capillary length k^^, we can approximate (/'(k) by (/'(O) = 1. In this limit, equation 
([l9| ) takes the following form: 

1 1 fV -Vr\ 

■\n{M~ 1) ^ — ^ha — . (20) 



27r%/27 27^^/27 V Vc 

If we combine equations (|l^) and (|20|), we see that slightly above the threshold, the wave- 
resistance behaves like: 



7f Vln' {V/Vc - 1) 



R--^[ ._.s:: .J - (21) 



Equation (^1|) constitutes the main result of this paper. First, we notice that for small 
objects this relation is, as expected, independent of the actual shape of the pressure field. 
Second, and more important, it shows the wave resistance R cancels out at y = Vc- This 
smearing is due to the cancellation of the vertical force Fq near the threshold that we get from 
the behaviour of S. 

Conclusion. - In this paper we have shown that in a fixed-depth situation the discon- 
tinuity of the drag-force calculated in |^ vanishes and is replaced by a smooth variation, in 
accordance with the experimental results found in |7| . A quantitative comparison between the 
present analysis and the data of Ref. |Q is however more involved since experiments from 
were performed in narrow channel geometry [ p^ with objects of size comparable with the 
capillary length. To recover the scaling relation observed experimentally, it would also be 
necessary to take into account both the variations of the shape of with V and the reflections 
of the waves on the walls of the channel. The present study suggests nevertheless that to fully 
test relation (H), experiments need to be devised that would measure both R and Fq. 
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